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Introduction. - Recently, there has been considerable 
interest in the investigation of critical random matrix 
ensembles (CrRME) (see [2-0 and references therein). 
CrRME are described by N x N matrices M mn whose 
characteristic feature is the slow decrease of off-diagonal 
matrix elements [J] 



(1) 



Such ensembles were introduced to model the Anderson 
transition of electrons in a disordered potential. The 
transition occurs between localized and extended states. 
Let A(a) and ^> m (a) (labeled throughout the paper by 
Greek letters) denote the eigenvalues and eigenfunctions 
of M mn . Localization properties of eigenfunctions can be 
described by a set of multifractal dimensions D q defined 
by 
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where (. . .) is the average over some eigenvalue window 
and over random realizations of the matrix. For localized 
states all D q are equal to 0, while for states delocalized 
over the whole d-dimensional space D q = d; for CrRME, 
states are multifractal and D q are non-trivial functions 
of q. Statistical properties of the eigenvalues of CrRME 
can be described by the level compressibility \- It ls 
defined from the limiting behavior of the spectral number 
variance 



E( 2 >(L) = (n(L) 2 )-(n(L)) 2 
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where n(L) is the number of eigenvalues in an interval L 
(the spectrum is unfolded with mean level spacing 1, so 
that (n(L)) = L). For the Poisson statistics of indepen- 
dent random variables \ — L while for standard random 
matrix ensembles x — ^ typically for CrRME one has 
< x < !• 

Multifractal dimensions D q are related with eigenfunc- 
tions, while compressibility x is related with eigenvalues. 
Thus there is no obvious relation between them. Never- 
theless in [H it was argued that for d-dimensional systems 
one should have 2x + D2/d = 1. Later, it was understood 



that this relation is valid, in general, only in a weak mul- 
tifractality regime, i.e. at first order in the deviation from 
the usual random matrix limit D 2 = d, x = Prig] . 

In this Letter we argue that, for three different one- 
dimensional CrRME considered below as well as for cer- 
tain two- and three-dimensional systems, the following 
relation holds: 



X + D 1 /d=l. 
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Here D± is the information dimension, corresponding to 
the mean eigenfunction entropy averaged over the same 
window as in 
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D\ In N. (5) 



We are not aware of general analytical arguments in favor 
of the conjecture (0]), as the fractal dimensions are not di- 
rectly accessible for analytical calculations. Nevertheless, 
a perturbation series approach provides an analytical way 
to them. There exist two regimes of perturbation series: 
strong multifractality when D q is closed to the Poisson 
value, D q <C 1, and weak multifractality when D q is near 
the random matrix value, d — D q « 1 [8]. Note that in 
the two extreme cases of Poisson and usual random ma- 
trices, ([4]) is trivially verified. For all CrRME with d = 1 
considered below we checked analytically that at first or- 
der of the perturbation series the fractal dimensions have 
the universal form 



D q = 



1 i - qx for x < 1 
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in a certain range of values of q. Thus (0} is valid at 
leading order of perturbation series. To check relation (HJ 
for intermediate values we performed careful numerical 
computations of both Di and x- The main result is that 
with available numerical precision no contradiction with 
our conjecture has been observed. 

Critical power law random banded matrices. - The 
most investigated CrRME is the ensemble of critical 
power-law banded random matrices (PLBRM) [lrS], @- 
8] . This is the ensemble of N x N matrices (real symmet- 
ric for f3 — 1 and complex Hcrmitian for (3 — 2) whose 
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matrix elements are independent random Gaussian vari- 
ables with zero mean and variance (depending on a pa- 
rameter b) given by (\H nn \ 2 ) — /3 _1 and for m ^ n 



(7) 



In both perturbative regimes of large and small b, the 
fractal dimensions and the level compressibility have 
been calculated at first order 0, ||, and it is easy to 
check that in these regimes ^ is fulfilled for this model. 
The second order terms at small b for \ an d -^2 have 
been calculated in 0| and [H respectively but the result 
for Di is yet unknown. 

To check the conjecture (0| for intermediate values of b 
we perform numerical calculations of D\ and \ for a crit- 
ical PLBRM where to reduce boundary effects the term 
m — n in is replaced by N/wsm((rn — n)Tr/N) Q. The 
fractal dimension D\ is extracted from a fit of the mean 
entropy ([5]) of the form a+Di \nN + b/N. The mean and 
variance of the entropy are calculated for eigenvectors of 
PLBRM of size 7V = 2",8<n<13. Average is per- 
formed over 8192 eigenvectors (namely 2™ -3 eigenvectors 
with eigenvalues around the band center, and 2 16 ~ n real- 
izations of the random matrices). The number variance 
([3]) is calculated on the unfolded spectrum with mean 
level spacing A = 1 by taking the average over windows 
of length L = 2kA, 1 < k < 32, centered at integer 
positions of the energy E = — 32 to E = 32, and over 
r realizations of the random matrices (from r = 32000 
for N = 256 to r = 500 for n = 2048). The level com- 
pressibility x(N) is then extracted from a quadratic fit 
£( 2 )(X) = a + x(N)L + cL 2 in the range L G [10,32]. 
The large- N asymptotics for \ is finally obtained by a 
linear fit of x(iV) as a function of 1/N over the range 
2 8 < N < 2 n . 

The results for Hermitian matrices (/3 = 2) are pre- 
sented in Fig. [TJ The agreement between \ and \ — D\ is 
quite good for all b. We obtain similar results for the real 
symmetric case /3 = 1 (data not shown), which indicates 
that indeed ((U holds for PLBRM matrices. 

Ruijsenaars-Schneider ensemble. - Our second ex- 
ample of CrRME is the Ruijsenaars-Schneider ensemble 
(RSE) proposed in [lOj. This ensemble consists of unitary 
matrices related to the Lax matrix of the Ruijsenaars- 
Schneider classical iV-body integrable model. Matrices 
of this ensemble have the form 
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1 < m, n < n. Here $ m are independent random vari- 
ables (phases) uniformly distributed between and 27r, 
a is a free parameter independent on N. The advantage 
of this model is that all spectral correlation functions, in 
particular the spectral compressibility \> can be calcu- 
lated analytically for N -> oo [IFJ, [ll|. When < a < 1, 




FIG. 1: (Color online) \ (black circles) and 1 — D\ (red error 
bars) for PLBRM Hermitian ensemble (,8 = 2). Straight lines: 
asymptotic theoretical values x = 1/ (47rfc) for 6 > 1 (solid) 
and x = 1 ~ Trv^fe - 4(2/^ - 1)tt 2 & 2 for b < 1 (dashed) 
(b 2 term is taken from pj and corrected for a misprint [3]). 
Circles for \ are larger than error bars. Inset: same data 
zoomed in. 



x(a) = (1 — a) 2 [11]. The general expression for \ = x( a ) 
is tedious and will not be presented here. We only men- 
tion that in the vicinity of an integer k 

(a-k) 2 

In order to obtain analytical expressions for D\ 1 
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we construct the perturbation series around all integer 
points a = k. Let a — k + e and expand the matrix 
elements ((8]) into series of e. One gets 



MM 1 + 
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where Mil = e'* m (5„, m+ fc, and 
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N sin(7r(m - n + k)/N) 
1 when n = m + k mod TV and oth- 



(here 8 n m +k — l wnen n = m 

erwise). In the case k = (i.e. \a\ <C 1), Mml is diag- 
onal: the unperturbed eigenfunctions are localized, and 
the perturbation decays as in (TT]). The first order of the 
perturbation series can be obtained by calculation of con- 
tributions from 2x2 sub-matrices @, H J8 1 , and one can 
check that in this case © is valid (see ll| for details). 

For other values of fc, eigenfunctions of the unper- 
turbed matrix Mmn are extended and the perturbation 
series corresponds to weak multifractality. Perturbation 
series in this regime have been constructed for critical 
PLBRM by using the supermatrix a- model 0, Q- For 
RSE, the form (flU]) of the matrix enables us to use 
standard perturbation series formulas. Let us consider 
the case k = 1. Eigenvalues A'°)(a) and eigenfunctions 



M?(a) of the unperturbed matrix are given by 
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where $ = J2f=i ®j/ N and 



5 n (a) = ^a(n-l)-^(^--$). (13) 
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The expansion of the exact eigenfunctions into a series of 
unperturbed eigenfunctions has the form 



N 

*„(a) = *( )( a ) + ^C Q ^(0) (/3 ). 
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At first order in e = a — 1, 
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A(°)(a)- A(°)(/3) 
Expansion at leading order in e yields 
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Using the explicit expressions (fT2"j) - (fT3")) one finds that the 
exact second order contribution to IF(cv) is 



JV-l N-l 



W(a ) = £ 2 — VV Sin in(3n/N) 

N3 t=\ sin 2 (vm/iV) sin 2 (tt/3/ A) 
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We are interested in its behavior for N —> oo. The only 
diverging terms correspond to two regions. The first is 
/3 -C N with n/N of the order of 1 and the second is 
n -C N and f3/N ~ 1. In this approximation 
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with y = n/N and < y < 1. The remaining sum over n 
can be transformed into an integral over variable y and 
finally we obtain 



W{a) - 2e 2 ln7V + C(l). 

N— >oo 



(21) 



From p^|) and © it follows that in the leading order of 
perturbation series in e = 1 — a one has 



A, = 1 - <?(1 - a) 2 



(22) 




FIG. 2: (Color online) \ (solid line) and 1 — D\ (red circles) 
for RSE. D\ is obtained numerically by averaging over all 
eigenvectors taken from 128 realizations for N — 2 s till 8 for 
TV = 2 12 . For X we use the theoretical value obtained in 
[r0 |. Dashed line: perturbative regime Inset: zoom for 
1< a < 3. 



For k > 2 calculations are more tedious but one can show 
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that when la — A;| <C 1 



D q = l-q 



(a - kf 
k 2 
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Comparing with @ we conclude that the leading terms 
of perturbation series in RSE indeed verify © around 
all integer values of a. In Fig. [5] we show that x = 1 — 
D\ is fulfilled for other values of a as well, with good 
precision. A variant of RSE has been investigated in 
fl2j ] and the same relation ^ has been observed to hold 
within numerical precision. Note that the perturbation 
series approach used here for weak multifractality can be 
applied to other problems as well (Tlj . 

Critical ultrametric ensemble. — As a third exam- 
ple we consider the ensemble of critical ultrametric ran- 
dom matrices proposed in [l3j]. This ensemble consists 
of 2 K x 2 K Hermitian matrices whose matrix elements 
are independent Gaussian random variables with zero 
mean. All diagonal elements have the same variance 
(\H nn \ 2 ) — W 2 . The variances of off-diagonal elements 
are (\H mn \ 2 ) = 2 2 ~ dmn J 2 , where d mn is the ultrametric 
distance between m and n on the binary tree with K lev- 
els and the root at 1. The parameter in this model is the 
ratio J/W. The first term of perturbation series in J/W 
for fractal dimensions has been calculated in 13j as 
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The calculation of the spectral compressibility can be 
performed similarly as above llj and one finds 



x = i- 
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Thus the relation x — 1 — -Di is fulfilled for the critical 
ultrametric ensemble at first order in J/W . For other 
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FIG. 3: (Color online) Same as Fig.[T]for ultrametric matrices. 
Solid line is the asymptotic theoretical value (|25[) . Dashed 
horizontal line is the numerical value 1 — D\ ~ 0.137 ± 0.008 
obtained for W — (J/W — > oo). Same method as in Fig. Q] 
(here Di is obtained from matrices of size up to 2 12 only). 
Inset: same data magnified around J/W = 0.5. 



J/W, we have calculated D\ and x numerically. The 
results shown in Fig. [3] confirm that the relation (Q| is 
valid with a good accuracy for all values of J/W. 

Higher dimensional models. - The above discus- 
sion was restricted to one-dimensional systems. Let us 
now turn to two examples of higher-dimensional systems 
where numerical data are available in the literature. 

The standard two-dimensional critical model is the 
metal-insulator transition in the quantum Hall effect, 
modeled by the Chalker-Coddington network [LJ]. We 



get the information dimension D± from 15|, where mul- 
tifractal dimensions for this model have been carefully 
fitted. Assuming that errors in the coefficients of the fit 
are independent we obtain that Di = 1.7405 ± 0.0004. 
According to our conjecture, the compressibility should 
be x c = 1- A/2, which givesy c = 0.1298±0.0002. This 
agrees with the estimate in [16[ where for this model the 
value x = 0.124±0.006 was obtained numerically. Notice 
that errors are underestimated as they mostly take into 
account statistical errors. 

The metal-insulator transition in three-dimensional 
Anderson model is the most important example of criti- 
cal systems. In [TtJ it was reported that Di = 1.93±0.01. 
Assuming that the s ym metry for the multifractal spec- 
holds, then the estimate D' n = 



trum conjectured in 
4.027 ± 0.016 given in 



[18| yields D 1 = 2d - D' a 



1.973±0.016. According to our conjecture @ one should 
have Xc = 1 — -Di/3, which gives Xc ~ 0.34 to 0.36. 
The spectral compressibility for the anisotropic Anderson 
model at the metal-insulator transition has been reported 
in [2(J. In this paper it was concluded that x = 0.28±0.06 
but this value corresponds to the average over different 
fits which have big fluctuations. In the same paper, when 
a smaller length has been used to define the number vari- 
ance it was found that x fluctuates much less and gives 
X = 0.32 ± 0.03, which agrees with Xc ~ 0.35 obtained 
from our conjecture. It would be of interest to get D\ 



and x with higher precision for these models. 

Conclusion. - In this Letter we present analytical and 
numerical evidences in favor of the conjecture that for a 
large class of CrRME the wavefunction entropy and the 
level compressibility are simply related by x + D\/d = 1. 
We consider three different models: the standard crit- 
ical PLBRM the ensemble of random matrices re- 
lated with Lax matrices of the Ruijsenaars- Schneider in- 
tegrable model [HI], an d the critical ultrametric ensem- 
ble [HI]. For all these models we check the conjecture in 
perturbation series and by direct numerical calculations. 
Though we cannot rigorously prove our relation, these 
investigations show that with good numerical precision 
it is fulfilled for very different systems. This suggests the 
existence of an universal structure in CrRME. 
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